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ABSTRACT 

We observed the oxygen-rich Large Magellanic Cloud (LMC) supernova remnant N132D (SNR 0525- 
69.6), using all instruments onboard the Spitzer Space Telescope, IRS, IRAC, and MIPS (Infrared 
Spectrograph, Infrared Array Camera, Multiband Imaging Photometer for Spitzer) . The 5-40 /im IRS 
spectra toward the southeastern shell of the remnant show a steeply rising continuum with [Ne III] 
and [O IV] as well as PAH emission. We also present the spectrum of a fast moving ejecta knot, 
previously detected at optical wavelengths, which is dominated by strong [Ne III] and [O IV] emission 
lines. We interpret the continuum as thermal emission from swept-up, shock-heated dust grains in the 
expanding shell of N132D, which is clearly visible in the MIPS 24 /im image. A 15-20 /mi emission 
hump appears superposed on the dust continuum, and we attribute this to PAH C-C-C bending 
modes. We also detect the well-known 11.3 /im PAH C-H bending feature, and find the integrated 
strength of the 15-20 /im hump about a factor of seven stronger than the 11.3 /im band in the shell 
of the remnant. IRAC 3-9 /xm images do not show clear evidence of large-scale, shell-like emission 
from the remnant, partly due to confusion with the ambient ISM material. However, we identified 
several knots of shocked interstellar gas based on their distinct infrared colors. We discuss the bright 
infrared continuum and the polycyclic aromatic hydrocarbon features with respect to dust processing 
in young supernova remnants. 

Subject headings: infrared: ISM — supernova remnants — supernovae: individual (N132D) 



1. INTRODUCTION 

Supernovae influence the chemistry and physics in the 
interstellar medium (ISM) from galactic scales down to 
the atomic level. Primarily, it is the energy input of 
supernova explosions that drives important processes 
in the ISM: strong shock waves accelerate, heat, and 
compress interstellar gas, destroy dust grains, acceler- 
ate cosmic rays, and affect interstellar chemistry (cf. 
IDraine fc McKellT^ . 

The role that supernovae play in the formation, de- 
struction, and processing of polycyclic aromatic hydro- 
carbons (PAHs) is largely unexplored. PAHs are believed 
to be the carrier of the unidentified infrared (UIR) bands, 
well-studied emission features near 3.3, 6.2, 7.7, 8.6, and 
11.3 /xm attributed to the C-C and C-H stretching and 
bend i ng modes of these m olecules l|Allamandola et aD 
119891 iPueet k. Legeil 11989(1 . It is established that in- 
terstellar gr ains are read i ly de stroyed by strong shock 
waves (e.g. iTielens et alJ ll994). and it was suggested 
that grain destr uction through sho cks may be a source of 
PAH molecules l|Jones et alJll996(l . Although PAH fea- 
tures have been observed in a wide variety of sources (cf. 
iVan Diedenhoven et alJl2004(l . no PAH emission has yet 
been detected in S NRs as far as we are aware, however, 
iReach et alJ l|2006|l identified four SNRs whose infrared 
colors suggest emission from PAHs. 

More recently, spectral features at 15-20 /im including 
broad plateaus and distinct, narrower features, most no- 
tably near 16.4 and 17.4 /tm, were detect ed in conjunc- 



and see iPeeters et alJ 12004 for a summary). The 15- 
20 /im features are generally attribut ed to C-C-C bend- 
ing modes of large PAH molecules ( cfJAJkimjmdola^^ilJ 
fT989t IVan Kerckhoven et alJl2000t IPeeters et all 12004(1 . 
The relatively narrow observed widths of the discrete 
16.4 and 17.4 /im bands suggest that they arise from in- 
divid ual PAH molecules in the gas phase ijPeeters et alJ 
12004(1 . There is also evidence that the appearance, rel- 
ative strengths, and profiles of the PAH features are 
coupled to the local physical conditions lIHonv et alJ 



coupled to tne loc al p nysical conditions (riony et^Lj 
20011: IPeeters et alJl200l IVan Diedenhoven et alJ 120041: 



Werner et ai] l2004). 
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In this paper, we report the detection of weak PAH fea- 
tures near 6.2 and 7.7-8.6 /tm, at 11.3 /im, and a promi- 
nent, broad hump at 15-20 /tm from the young supernova 
remnant N132D (SNR 0525-69.6) in the Large Magel- 
lanic Cloud (LMC). The age of N13 2D determined from 
the ejecta kinematics is ~2500yr (M orse et a[]|l995), 
and its spatial extent is about 100", or 25 pc assum- 
ing a LMC distance of 50kpc. N132D belongs to the 
small class of young, oxygen-rich SNRs that are pre- 
sumably the product of core-collapse supernovae, with 
the most famous Galacti c example b eing CasA. For 
the progenitor of N132D, IBlair et al . ( 2000) suggested 
a W/O star of 30-35 M Q based on relative elemen- 
tal abundances and comparisons to stellar nucleosyn- 
thesis models, which may have exploded as a Type 
lb su pernova. N132D show s a radio spectral index of 
-0.7 ((Dickel k Milnd lT9951. consistent with its overall 
shell-like morphology. It is one of the brightest X-ray 
sources in the LMC (IHughesI 119871: [Hwang et alJll993t 
iFavata et alll997llHughes et al.ll998(l . and h as also been 
investi gated in great detail i n the UV/ opti cal ( Blai r et all 
Il99i IMorse et all Il995l ISutherland fc Dopftal 1 1993 
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iMorse et all Il996t iBlair et all [20001. The large-scale 
morphology of the gas in N132D s hows three types o f 
line-emitting regions in the optical ijMorse et alj [l9961: 
(1) red- and blueshifted oxygen-rich ejecta filaments in 
the center, with a total velocity range ~ 4400 km s -1 
l)Morse et al.ll9 951. (2) shocked interstellar clouds within 
the remnant perimeter and in the outer rim, and (3) a 
photoionized precursor around the outer edge of the rem- 
nant. 

We observed N132D in the infrare d with the Spitzer 
Space Telescope ijWerner et al J 120041. In Section |2 we 
describe our observations and the data reduction proce- 
dure. We took particular care of the background sub- 
traction in order to avoid confusion with foreground and 
background emission ( Section 12. 2|) . We present our re- 
sults from the imaging and the spectroscopy in Sectional 
and give a detailed discussion of the IR emission features 
and a comparison to PAH features observed in other as- 
tronomical objects in Section 0] 

2. OBSERVATIONS AND DATA ANALYSIS 

We observed N132D with all instruments onboard the 
Syitzer Syac e Telescope: IRS (Infrared Spectrograph, 
iHouck et all 120041 2004 Decembe r 13, IRAC (Infrared 
Array Camera. iFazio et all l200l) 2004 November 28, 
and MIPS (Multi band Imaging Photometer for Spitzer, 
iRieke et all I2004H 2004 November 7 (Spitzer aorkeys 
11019776, 11018496, and 11022336, respectively). We 
used the IRS in its low resolution setting, yielding a 
wavelength coverage of 5.2-8.7/7.4-14.5 /im with the 
short-low (SL2/SL1) and 14.0-21.3/19.5-38.0 /xm with 
the long-low (LL2/LL1) modules. The nominal spec- 
tral resolution is R = 64-128 as a function of wavelength 
for each module and order, e.g. R ~ 100 at 18 fim. The 
spatial sampling of the IRS is 5.1/1.8" per pixel in the 
LL/SL modules, and the spatial resolution is diffraction 
limited by the 85 cm primary mirror. The four IRAC 
bands cover the wavelength ranges 3.2-4.0, 4.0-5.0, 5.0- 
6.4, and 6.4-9.4 /im, respectively, with a pixel size of 1.2". 
The three MIPS detector bands cover 20.8-26.1, 61-80, 
and 140-174 /im, with respective pixel sizes of 2.6, 5.2, 
and 17" in our configuration. The approximate total ex- 
posure times are 270 s for each IRAC band, 48, 38, and 
25 s for the MIPS bands, 280 s for each IRS SL, and 360 s 
for each LL module. 

2.1. General data reduction procedure 

The pipeline processed MIPS and IRAC image mosaics 
were of good quality, and we used them without further 
processing. The analysis of the IRS spectra is compli- 
cated by the fact that N132D is a non-uniform, extended 
source superposed with a complex background adjacent 
to a molecular cloud (see Fig.^l. Moreover, the Spitzer 
science pipeline is optimized for point sources, and the 
extraction of extended sources from the default pipeline 
products results in biased spectra. Our principal data 
reduction strategy for the IRS spectra consists therefore 
of the following steps: 

(1) starting from the S13.2.0 pipeline products, we in- 
terpolated bad pixels flagged by the pipeline, and re- 
moved rogue pixels by subtracting nearby sky positions 
taken from the Spitzer archive for each SL/LL module 
(see Sect. l2~2l for details). Any remaining rogues were 




Fig. 1. — N132D: Spitzer IRAC 3-9 /mi color composite image 
(3.6 /mi=blue, 4.5 /mi=green, 5.7/mi=orange, and 8.0/mi=red); 
contours at the image center are X-ray ROSAT counts of N132D, 
and contours south of the remnant ar e CO J = 2 — > 1 emission 
(both adapted from Banas et al. 1997). The IRS slit coverage is 
outlined with polygons: SL=west-east, LL=north-south; yellow 
(SL) and magenta (LL) boxes indicate the extracted parts for the 
southeastern rim spectrum, the green box marks the position of an 
ejecta knot detected in the IRS LL spectra, and blue boxe s mark 
the extracted positions for the local background (Sect, IT^l . 

removed with the IRSCLEAN 1 tool. Rogue pixels are 
single detector pixels that show time-variable and ab- 
normally high flux values. 

(2) we extracted the spectra with SPICE 2 v.l.3betal 
using the default extr action widths, and imported the 
spectra into SMART 3 l|Higdon et al.ll2004H for postpro- 
cessing; 

(3) subtraction of a local background spectrum (see 
Sect. 152)1: 

(4) basic aperture and slit loss correction assuming a 
uniform, extended source, to account for the bias in- 
troduced by the point source optimized Spitzer science 
pipeline. We derived a correction function by modelling 
the Spitzer point-spread-function (PSF) for a range of 
wavelengths with the software STINYTIM 1 developed by 
John Krist, and then dividing the fraction of the PSF ad- 
mitted by the slit as defined by the slit-width and wave- 
length dependent extraction aperture by the sum over 
the whole PSF. The resulting flux correction functions 
have wavelength dependent values of 0.8-0.6 for SL and 
0.85-0.65 for LL. Further division of the corrected flux 

1 available through the Spitzer Science Center at Caltech, 
see http : //ssc . spitzer . caltech . edu/archanaly/contributed/ 
IRSC-LKAJN is developed by Jim ingalls with contributions from 
the Cornell IRS instrument team. 

2 Spitzer IRS Custom Extraction software available 
through the Spitzer Science Center at Caltech, see 
http : //ssc . spitzer . caltech. edu/postbed/ spice .html 

J The Spectroscopy Modelling Analysis and Reduction Tool was 
developed by the IRS Team at Cornell University, and is available 
through the Spitzer Science Center at Caltech. 
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density by the extracted area given by the slit-width and 
the extraction aperture yields spectra in surface bright- 
ness units, i.e. flux density per steradian, which can be 
directly compared to the IRAC and MIPS fluxes. Note 
that IRAC fluxes are calibrated for point sources, but es- 
tablished point -to-extended sour ce flux correction factors 
are available ijReach et alJl2005^ : 

(5) we removed fringes that appear in LL1 at 19.5- 
24.0 /im with the SMART defringing tool, averaged mul- 
tiple spectra of the same region, and finally merged SL 
and LL spectra. Note that we extracted the SL and LL 
spectra of the southeastern rim from two slightly offset 
positions, because the small overlap region covered by 
both slits is positioned just outside the brightest regions 
in the southeastern rim (cf . Fig. . The two extracted 
regions are separated by about 10", and the MIPS 24 /mi 
fluxes in both regions differ by only 10 percent, whereas 
the IRAC channel 3 and 4 fluxes are comparable. 

2.2. Background removal 

A meaningful background subtraction of the spectra is 
one of the most critical and difficult tasks, and of ma- 
jor importance for this work. Commonly, either a spa- 
tially close spectral nod or a dedicated off-source position 
are subtracted to remove the background. Both of these 
methods are problematic in our case, since we are facing 
a complex, non- uniform background (cf. Fig. ^ . 

As a first step, we subtracted archival off-source sky 
exposures from all spectral modules (cf. Sect. 12.11 item 
1). The selected sky positions are devoid of continuum 
and line emission, and were observed close in time and 
position with respect to our observations (< 3 hrs, ~ 10'). 
The sky position Spitzer aorkeys are 10964224, 12935424, 
and 12936960 for SL2, 12936960 for SL1, and 12935424 
for LL2/LL1. The sky subtraction removes rogue pix- 
els, zodiacal light emission, and flat-field imperfections, 
which are particularly noticeable in the SL modules. 

In a subsequent step, we subtracted a rogue-cleaned lo- 
cal background spectrum from the extracted source spec- 
tra for each spectral module (cf. Sect. 12. ll item 3). We 
selected positions that lie close enough to the source to 
represent the local background but sufficiently far from 
the rim to avoid the influence of the radiative precursor 
that extends noticeably beyond the edge of the remnant 
(see Sect. l3.l|) . The background positions are depicted as 
blue boxes in Fig.^ We averaged two background posi- 
tions from opposite sides of N132D for SL2 and LL2/LL1, 
while ensuring that both positions show similar spectra. 
For SL1, we only have slit coverage east of the remnant, 
and we adopted the same background position as for SL2. 
The overall background spectrum extracted from these 
positions is shown in Fig. 0] We estimate a 30 percent 
uncertainty in this background spectrum due to spatial 
variations by using slightly different extraction positions 
(approximately ±10"; blue boxes in Fig.QJ and two back- 
ground positions on opposite sides of the remnant. 

2.3. Analysis of emission lines 

We measured emission line fluxes, central wavelengths, 
and full-widths at half-maximum (FWHM) from line fit- 
ting with the IRAF 'spectool' package 4 . Due to the 

4 IRAF (Image Reduction and Analysis Facility) and 'spectool' 
are maintained by the IRAF programming group at the National 



low resolution of the IRS spectrograph in our setting, 
B ~ 100 at 18 /mi, the lines are usually unresolved 
or marginally resolved, and represent a convolution of 
the true line profile with the instrumental profile. We 
used Gaussians to fit the line profiles of atomic lines/H2 
and Lorentzians for discrete PAH features. Our quoted 
FWHM are not deconvolved from the instrumental IRS 
profile. 

We also extracted spectra from different spatial posi- 
tions at « 25" intervals along the LL spectrograph slit in 
order to study the spatial variation of line fluxes across 
the supernova remnant and the molecular cloud to the 
south (see Sect. ESI Fig.©. For this particular task, we 
employed a slightly different data extraction strategy. In 
order to improve the spatial resolution, we narrowed the 
default SPICE extraction aperture by a factor of two, 
giving extraction widths of about 12" for LL2 and 20" 
for LL1 near the central wavelengths of each module. 
This procedure systematically reduces measured fluxes 
by about a factor of two compared to the calibrated de- 
fault extractions, which is of no concern as long as we 
compare only relative fluxes of extended emission ex- 
tracted in that same fashion. We subtracted the con- 
tinuum by fitting a low-order Legendre polynomial using 
spectral regions that are free of emission lines, usually 
below 15 and above 19 /mi. For the southeastern rim 
of N132D, this method gives a similar continuum fit for 
the LL2 14-21 /mi region than the modified-blackbody 
fit shown in Fig. We did not apply a local background 
subtraction here, which means that we compare emis- 
sion line fluxes including fore- and background for each 
spatial position along the LL slit. 

3. RESULTS 

3.1. IRAC and MIBS imaging 

N132D shows strong emission in the MIPS 24 /im band 
with a background subtracted peak surface brightness of 
50MJy/sr _1 in the southeastern rim. The total, back- 
ground subtracted MIPS 24 /im flux integrated over the 
whole remnant is 3.0 ±0.2 Jy. The surface brightness dis- 
tribution follows the X-ray contours very well (cf. Fig. 
left panel), although note that the MIPS 24 /im emis- 
sion visibly extends beyond the edge of the remnant as 
defined by the X-ray contours, particularly along the 
western, southern, and eastern edges. This is a clear 
sign of the radiati ve precursor also seen in the optical 
ijMorse et al.l ll996), and is most likely continuum emis- 
sion from heated dust and/or line emission from ionized 
oxygen caused by hard photons from the remnant; it can- 
not solely be explained by the wider MIPS point-spread- 
function (cf. point source in Fig.|21 left panel). The rem- 
nant is also marginally detected in the MIPS 70 /im band 
with about 60 MJy/sr -1 background subtracted peak sur- 
face brightness and a total integrated flux of 9.6 ± 1.7 Jy, 
but the background confusion is much stronger (Fig. 
right panel). There is no conclusive detection in the 
MIPS 160 /im band, with an upper limit of 5MJy/sr _1 
peak surface brightness. The four IRAC bands do not 
show clear evidence of emission spatially correlated with 
the remnant as seen in MIPS nor with the ejecta seen in 
the optical (cf. Fig. |3 left panel) . This is partly due to 

Optical Astronomy Observatory (NOAO) in Tucson, Arizona. 



4 



Tappe, Rho, & Reach 




(MI/,'h1 twyxt 

Fig. 2. — N132D: (left panel) Spitzer MIPS 24 /an image superposed with archival Chandra ACIS 0.3-10 keV contours; note the star 
at 05 h 25 m 16?7, -69°38'39'/9 (classified as a semi-regular pulsating star bv lCioni et"ani200!l but being a source of hard X-ray emission), 
which shows the approximate width of the MIPS point-spread-function (PSF). The pixel scale is 2.55"/pixel or 0.62pc/pixel; (right panel) 
Spitzer MIPS 70 fim image superposed with Chandra ACIS 0.3— lOkeV contours. The pixel scale is 5"/pfxel or 1.21 pc/pixel. 




Fig. 3. — N132D IRAC emission knots: (left panel) Spitzer IRAC 3.2— 6.4 (im smoothed color composite (3.6 /im=blue, 4.5 ^tm=green, 
5.7/im=orange). The dashed line represents the outermost ROSAT X-ray contour from Fig. a nd the polygons represent the IRS slit 
coverage. The dotted contours outline the blue and redshifted ejecta of N132D as seen in the optical I Morse ct al. 1995); (right panel) Spitzer 
IRAC smoothed color difference image: 3.6 /xm=blue/ white (faint/bright stars), 4.5 fira— (3.6 /xmx0.65)=green, 5.7/^m— (8.0/^mx0.33) — 
(3.6 /xmx0.5)=red. The contrast levels have been c hosen to suppress background noise fluctuations such that only the strongest emission 
peaks appear against a black background. See Sect. l4~T1 for a detailed discussion. 



confusion with the foreground and background ISM ma- 
terial. The IRAC emission within the remnant perimeter 
is clumpy and shows voids that may have been created 
by the passing Shockwave of N132D. The IRAC color ra- 
tios measured at different positions within the remnant 
are consistent with PAH emission Ccf. lReach et~aT1 l2006. 
Fig. 2). We also detected several small emission knots 
that show IRAC colors different from the surroundings. 
In order to remove background stars and confusing PAH 
emission, we subtracted a scaled IRAC 3.6 ^m image 
from the 4.5 /im image, and a scaled IRAC 3.6 and 8.0 /im 
image from the 5.6 ^m image (cf. Fig. |3 right panel). 
We determined the optimal scaling factors empirically to 
produce a generally flat background. The emission knots 
have bee n observed previou s ly in the optical with HST 
WFPC2 l|Morse et al.lll995l 11996]) . and were identified 
as clumps of shocked, interstellar/circumstellar material 
(cf. Sect. I4.fl for a more detailed discussion). Generally, 



the strongest emission in all four IRAC bands is associ- 
ated with the molecular cloud/H II-region located south 
of the remnant (cf. CO contours in Fig. QJ. 

3.2. IRS spectroscopy 

Figure 0] compares the IRS spectrum extracted from 
the southeastern rim of NI32D with the local background 
spectrum. Figure shows the background subtracted 
IRS spectrum from the southeastern rim of N132D to- 
gether with a two-component modified-blackbody fit. 
Given the difficult background removal and the slightly 
offset extractions of SL and LL spectra (cf. Sect. I2.2|l . it 
is worthwhile to note that the spectra from all modules 
seamlessly connect, and the fluxes are consistent with 
background subtracted flux estimates based on the cali- 
brated IRAC and MIPS data. 

The modified-blackbody fit to the continuum uses 
an emissivity for 0.1 /im silicate/graphite particles with 
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Wavelength Cum) 

Fig. 4. — N132D southeastern rim and local background Spitzer 
IRS spectrum (online color version: magenta=SL2, red=SLl, 
green=LL2, blue=LLl; blue dashed=background) : The upper 
spectrum is without, and the lower spectrum is with the local back- 
ground subtracted (solid black lines). The dashed line is the local 
background spectrum extracted from the blue boxes in Fig.0 The 
SL parts of the spectra have been slightly smoothed for illustration 
purposes. 




5 10 20 50 

Wovelength (^m) 



Fig. 5. — N132D southeastern rim Spitzer IRS spectrum (online 
color version: magenta=SL2, red=SLl, green=LL2, blue=LLl): 
diamond data points with error bars and bandwidths are back- 
ground subtracted IRAC 5.7, 8/^m and MIPS 24, 70 /im fluxes 
measured at the position of the extracted IRS LL spectrum. The 
SL parts have been slightly smoothed (see Fig. |5] for the origi- 
nal, unsmoothed version). The solid and dashed lines represent a 
two-component modified-blackbody fit to the spectrum; the fitted 
temperatures are 58 and 110 K (see text for details). 

a silicate-to-graphite ratio of 0.2, following the ISM 
silicon and carbon elemental abundances derived for 
N132D by HI air et alJ l)2000|) and silicon abund ances 
derived from X-ray data bv iHughes et al.l l)1998j) and 
iFavata et alJ l|1998|) . The shape of this simple emis- 
sivity function is very similar t o the emissivity from 
the ge neral LMC dust model of iWeingartner fc Draind 
( 2001) in the 5-35 /nm range. Note the shoulder near 
10 ^tm in the modified-blackbody fit due to the sili- 
cate part of the emissivity, which reproduces the con- 
tinuum flux of our spectrum and therefore provides di- 
rect evidence that the heated dust in N132D contains 
at least some fraction of silicates. Modelling of the 
UV to near-IR extinction curve for LMC sightlines l)Peil 




5 10 15 20 25 30 

Wovelength (^m) 



Fig. 6. — N132D southeastern rim spectrum (color online version: 
magenta=SL2, red=SLl, green=LL2, blue=LLl): the local back- 
ground and the modified-blackbody continuum from Fig. [5] have 
been subtracted. The total error includes the reduction pipeline 
delivered errors and an estimate of the variation in the local back- 
ground (see text for details). 

19921 iMaccioni fc Perinottol 119941: iClavton et all l2003t 

Cartledge et al 2005) generally seems to indicate a larger 
silicate content in dust grains, i.e. favouring larger sili- 
cate to graphite abundance ratios. For comparison, we 
also used the emissivity from the LMC dust model by 
IWeingartner fc Draind l)200lD and alternative silicate-to- 
graphite abundance ratios of 1 and 4 for our 0.1 /zm sil- 
icate/graphite dust emissivity. The original emissivity 
data for silicate an d gr aphite particles are ta ken from 
iDraine fc Led" lfl98l and lLaor fc Draind ifTiffl! 5 . We de- 
liberately excluded the PAH component from the emis- 
sivity in order to analyze these features later in the 
continuum subtracted spectrum (Fig. EJ). The SMART 
package uses a modified-blackbody fitting function of the 
form 

Flux density = (l - c at ^ x ) ^ n t B x (T t ), (1) 

i 

where aQ\ is the normalized emissivity Q scaled by the 
factor a, and HB\(T) is the Planck function at tempera- 
ture T scaled by the factor fi, which accounts for the size 
of the emissive region. With the fitted free parameters a, 
fi, and T, all emissivity functions mentioned above pro- 
duce nearly equally good fits to the spectral data, and 
the fitted values for T and f2 are similar for all models. 
For a two-component fit with a silicate-to-graphite ratio 
of 0.2, we obtain temperatures of T\ = 58 ± 15 K and 
T 2 = 110 ± 6 K, and a ratio fii/fi 2 ~ 15. 

3.2.1. PAH features from the supernova blast wave 

Figure H3 shows the IRS spectrum of the southern 
rim. The foreground/background and the fitted dust 
continuum have been subtracted (cf. Sect. 12.21 and I3.2jl . 
The quoted total error consists of the pipeline-delivered 
statistical ramp uncertainties, combined with an esti- 
mate of the error due to the local background variation 
(cf. Sect. 12. 2|) . We detect several emission lines from 
atomic ions, tentative weak PAH features near 6.2 and 

5 Emissivity data files are available from the homepage of Bruce 
T. Draine, http://www.astro.princeton.edu/~draine 
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7.7-8.6 /im, a weak PAH 11.3 /jm feature, and a promi- 
nent, broad hump at 15-20 /im. The ionic emission lines 
are from Ne 2+ , 3+ , and S 2+ ions. We can not dis- 
tinguish between [O IV] 25.91 /im and [Fe II] 25.99 /im in 
our low-resolution spectra, but argue that most of the 
emission is from 3+ (cf. Sect. I3.3|l . The PAH features 
at 6.2, between 7.7 and 8.6 /im, and at 11.3 /im are the 
well-studied PA H C-H in-plane and out-of-plane bend- 
ing modes (e.g. IHonv et all 120011 IPeeters et all I2002t 
iVan Diedenhoven et alJl2004|) . The presence of discrete 
6.2 and 7.7-8.6 /im bands is uncertain due to the low 
quality of our spectra in that region (cf. total error 
in Fig. There is probably an underlying emission 
plateau at 6-9 /im, which has been observed in other 
sources and is probably associated with broad continua 
due to larg e PAHs, clusters of small PAHs, or very small 
grains (e.g.lAllamandola et alJl!989t IPeeters et a11l2002t 
iRapacioli et al.ll2005|) . ~ 

The broad hump at 15-20 /im is presumably due 
to PAH C-C-C in - and out-of-plane bending modes 
iAllamandola et all Il98l IVan Kerckhoven et all 1200(1 
IPeeters et all 12004 ). Figure shows the promi- 
nent emission hump between 15 and 20 /im in detail. 
The spectrum without background subtraction clearly 
shows the 16.4 /im PAH feat ure (cf. iMoutou et alJl2000l 
IVan Kerckhoven et alJ I2000I) . and additional features 
near 17.1 and 17.9 /mi. These distinct bands with widths 
~ 0.5 /im may arise from smaller PAHs in the gas phase. 
However, the individual bands are hardly distinguish- 
able after background subtraction, indicating that these 
features are primarily due to the background (cf. back- 
ground spectrum in Fig. |3J . 



TABLE 1 

N132D FIT DATA WITH ERRORS GIVEN IN PARENTHESES: 
SOUTHEASTERN RIM (FlG.|§J, LOCAL BACKGROUND (FlG.IH, AI\'I> 
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FlG. 7. — N132D southeastern rim spectrum (LL2): detailed 
view of 15—20 /im region from Fig. El F° r comparison, we show the 
spectrum without the local background subtraction (dot-dashed 
line), and the total error from Fig. |5]as a gray-shaded error band. 

Table ^ lists the fit parameters for all detected ionic 
lines and PAH features in the rim, the local background, 
and the molecular cloud/H II-region near N132D. For the 
local background spectrum, we assumed a fiat continuum 
of 1.2 x lO^Wm-^m-'sr 1 . The detected PAH fea- 
tures in the southeastern rim of N132D are the first evi- 
dence of PAH molecules surviving a supernova blast wave 
to our knowledge. PAH features also appear in lines-of- 
sight outside the remnant perimeter, and it is necessary 



to make certain that the PAH emission originates from 
the blast wave, and not from foreground, background, or 
ambient ISM material. There are three points of evidence 
that the PAH emission is indeed originating from mate- 
rial associated with the blast wave. Firstly, we removed 
the foreground and background emission by subtracting 
the flux from lines-of-sight just outside the blast wave 
perimeter (Sect. l2~2"|) . Secondly, we analyzed the spatial 
variation of emission line fluxes along the LL spectro- 
graph slit (cf. Sect. 13.31 Figure EJ, which shows that the 
PAH features between 15 and 20 /im have a well-defined 
local flux maximum at the position of the rim. This is 
expected for material that has been swept up by the blast 
wave and is not yet completely destroyed by the shock. 
If the majority of the flux would be from unrelated fore- 
ground or background emission, there would be no local 
maximum at the position of the rim. Some fraction of 
the detected PAHs may reside in the radiative precursor 
(cf. Sect. 13.1( 1. where they would be excited/processed 
by energetic photons emanating from the supernova rem- 
nant. Finally, the 15-20 /im PAH emission changes no- 
tably when comparing features from the supernova rem- 
nant to the adjacent molecular cloud/H II-region and 
the local background (cf. Fig. W\ and ITT|l . which we in- 
terpret as PAH/grain processing in N132D (cf. Sect. 14.31 
and EXT}. 

3.3. Spatial variation of the emission features 

In Figure [S] we show the spectrum of the southeast- 
ern rim of N132D in comparison to the ejecta knot de- 
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N132D SE-Rim 




Wovelength (^m) 

Fig. 8. — Comparison of different Spitzer IRS spectra for N132D, 
its nearby molecular cloud/H n-rcgion, and an archival spectrum 
of the galaxy NGC 7331 (dashed line). We extracted and processed 
the average spectr um of the nucleus and inner ring of NGC 7331 
ISmith et al.ll2004l) in the same way as our spectra, and app lied 
a sca ling factor of 0.5 for illustration purposes (see Sect. EOl and 
I4.3.1l for further discussion). 




Projected distance from shock (pc) 

Fig. 9. — Spatial variation of line vs. PAH 15—20 (im fluxes, and 
total integrated flux from continuum and lines in the 14—21 /im 
range as a function of projected distance from the supernova shock. 
[Sill] and [O iv] appear only in LL1 (19.5—38 /im), limiting the spa- 
tial coverage to about 20 pc. The H2 flux across the supernova 
remnant is uncertain due to strong blending with PAH features 
but most likely very small (see text for further discussion). 



tcctcd in the LL spectra, the molecular cloud/H II- 
region south of N132D, and the average spectrum of 
the galaxy NGC 7331. The characteristic features of the 
rim spectrum have been described earlier in this paper. 
The spectrum of the ejecta knot lacks the steeply ris- 
ing dust continuum and the 15-20 /im PAH features, 
but shows strong ionic emission lines from Ne 2+ and 
3+ , which is consistent with the strong neon and oxy- 
gen lines seen in the o ptical for the same knot (N132D- 
P3, IBlair et all 12000]) . The line fluxes and FWHM 
are 1.6 x 10 _8 Wm _2 sr _1 (0.15 /im) for the neon and 
8.4 x 10~ 9 Wm~ 2 sr _1 (0.34 /zm) for the oxygen line. 
The spectrum of the molecular cloud/H II-region shows 
strong PAH bands as well as ionic and H2 lines. The 
association of the CO emission south of N132D with 
an H I I-region seen in Ha was noticed by iMorse et al.l 
(1995). iBanas et"afl l)1997t) estimated the mass of the CO 
cloud to be about 3 x 10 5 M Q , and argued that its spa- 
tial proximity and matching LSR velocity provide strong 
evidence that this giant cloud is associated with N132D. 
This provides an opportunity to relate our observations 
of dust and PAHs in the supernova remnant N132D to 
the nearby environment its progenitor star probably orig- 
inated from. Such a comparison will also be useful to in- 
terpret spectroscop ic observations of s tar-formi ng gal ax- 
ies (e.g. NGC 7331. ISmith et al.ll200l see Sect. Oil . 

We present a systematic assessment of the PAH 15- 
20 /im, ionic, and molecular line fluxes from positions 
across N132D and the nearby molecular cloud in Fig- 
ure We extracted spectra along the LL slit with 
« 25" separation starting at the northern edge of N132D 
and ending at the southern edge of the molecular cloud 
about 3.3' south of N132D, resulting in a total number 
of 15 spectra. We converted the spatial distance along 
the slit into a parsec scale assuming a LMC distance of 
50 kpc, and adopted the location of the southeastern rim 
of N132D as the zero position. Each node in Figure 
represents a line flux derived from continuum and line 
fitting, with errors given for the positions covering the 
supernova remnant (cf. Sect. 12. 3|) . 



The ionic and molecular lines detected and exam- 
ined in our spectra are [Ne III] 15.6, [S III] 18.7 and 
33.5, [O IV] 25.9, [Si II] 34.8, and H 2 0-0 S(0) 28.2 and 
S(l) 17.0 /im. The lines from atomic ions, particularly 
the strong [Ne III] and [S III] transitions, are mostly un- 
blended and straightforward to fit, which makes them 
excellent references in terms of flux and central wave- 
lengths. 

[Ne III] and [O IV] lines only appear across the su- 
pernova remnant and in the precursor, and are particu- 
lary strong at the position of the ejecta knot. Generally, 
[Ne III] and [O IV] trace the hot ionized gas of N132D. 
Note that we can not distinguish between [O IV] 25.91 /im 
and [Fe II] 25.99 /im in our low-resolution spectra, but we 
argue that most of the flux is likely to be from [O IV] 
because the flux follows [Ne III] closely (cf. Fig. ; un- 
der the conditions where we observe Ne 2+ , iron would 
be mostly present as Fe 3+ and Fe 4+ . There is a weak 
line near 17.9 /im, which may be [Fe II] 17.9, but pos- 
sible confusion with the PAH 15-20 /im features renders 
this line inconclusive. We do not see [Fe II] 24.5/35.3 /im, 
[Fe II] /[Fe III] 22.9, and [Fe VI] 19.6 /im, which are prob- 
ably too weak to be detected in our spectra. However, 
X-ray emission from Fe 16+ and more high ly ionized st ates 
was detected in the southern rim (|Behar et alJl200l|) . 

[S III] is a common line seen in PDRs and H Il-regions. 
Here, it not only appears across the remnant, but also 
across the molecular cloud tracing the H II-region. The 
[S III] 1-0 transition at 33.5 /im follows the same gen- 
eral trend, with the 2-1/1-0 ratio increasing from 0.08 to 
about 0.2 from the northern edge of the molecular cloud 
into the remnant. 

H2 traces the molecular cloud seen in CO (Fig.^J. We 
checked for H 2 0-0 S(0) 28.2, S(2) 12.3, and S(3)9.7//m, 
and did not detect these lines across the remnant; weak 
S(l) 17.0 /im emission is not detectable because of strong 
blending with PAH emission. We conclude that H 2 is 
mostly absent across the remnant, where it was prob- 
ably destroyed by the fast shock and the strong radia- 
tion field, but apparently present in the precursor. The 
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H2 S(1)/S(0) ratio increases from about 1.4 at the north- 
ern edge of the molecular cloud/H II-region to >11 at 
the shock precursor of the rim, which implies an increase 
in the UV radiation field as expected for the presence of 
a radiative precursor (cf. Sect. 13 . l|> . 

[Si II] shows a nearly constant flux, except a noticeable 
minimum at the position of the rim. This is most likely 
an effect of the increased radiation field and the elevated 
temperature. Silicon will be much more highly ionized 
under those conditions, which is confirmed by the detec- 
tion of X-ray emiss ion from Si 12+ in the southern rim by 
iBehar et al.l IpOOl). Unfortunately, this conceals a pos- 
sibly increased amount of silicon in the gas phase due to 
the destruction of dust grains. 

There are three well-defined PAH components at 16.5, 
17.4, and 17.8 /xm in the region of the molecular cloud 
(Fig. 181 and 1111 and Table 0, which have all been pre- 
viously observed in other sources (see Sect. 14 . 3TTfl . The 
fluxes of these features follow the H2 flux very well as ex- 
pected if the PAHs are well-mixed with the gas (Fig.EJ. 
The feature at 16.8 /xm is blended with the stronger PAH 
feature at 16.5 /im and the strong H2 S(l) line at 17.1 /xm, 
and therefore more uncertain. The situation changes 
for the supernova remnant, where the discrete 15-20 /xm 
PAH features all but disappear, and a nearly featureless 
15-20 /xm emission hump remains (Fig.0and Table ITjl. 

4. DISCUSSION 
4.1. IRAC emission knots 

We detect shocked interstellar clouds with IRAC due to 
their distinct emission in the IRAC channels (cf. Fig.[3J). 
The IRAC emission knots seen in the northern central 
part and the northwestern rim of the remnant coincide 
with previously detected, bright optical [O III] and [S II] 
emission complexes named "Lasker's Bowl" and "West 
Complex". They were identified as shocked interstellar 
clouds/circumstellar material based on optical line diag- 
nostics, elemental abundances, and their ra dial velocities 
similar to the local ISM surrou nding N132D ijMorse et al.l 
1199.1 119961: iBlair et all 12000ft . The small knot visible 
close to the center of the remnant is probably a smaller 
analogue of these two complexes. We do not detect the 
red- and blueshifted, oxygen-rich ejecta of N132D with 
IRAC (cf. the dotted contours in Fig. EJ> . 

The measured IRAC band ratios for the shocked clouds 
are roughly consistent with thermal emission from hot 
silicate/graphite dust grains. In particular, the West 
Complex shows an IRAC 4.5/8 /xm ratio of 0.18 and a 
3.6/5. 7 /xm ratio of 0.38. For this region, Bla ir et ahl 
( 2000) presented Hubble FOS spectra that show a broad 
range of ionization stages as expected from shock heat- 
ing, e.g. O to 3+ and N to N 4+ , and elemental abun- 
dances consistent with LMC interstellar material. They 
derive a high electron temperature T w 43, 000 K, a 
mean electron density of several thousand per cm -3 , and 
a preshock cloud density larger than 100 cm~ 3 . Their 
shock models indicate that a shock speed of at least 
~200kms _1 is needed to produce the observed high- 
ionization species. 

Comparing this to shocked d ust IR emission models 
calculated bv lDwek et alJ <|1996i). we find that the mea- 
sured IRAC color ratios are roughly reproduced by mod- 
els with a preshock density of 100 cm~ 3 and a shock 



speed <400kms _1 , and models with a preshock den- 
sity of 10 cm~ 3 and a range of shock speeds between 400 
and lOOOkms -1 . This suggests that part of the IRAC 
emission in the West Complex may be due to a popu- 
lation of shock-heated, small dust grains. However, it 
is likely that the measured IRAC colors represent a mix 
of dust continuum and ionic emission. The green chan- 
nel 2 excess emission probably includes H I Bret 4.05 /xm, 
since t he knots appear in the Ha image from lMorse et all 
( 1995) and the optical spectrum of the We st Complex 
show s strong hydrogen recombination lines ijBlair et alJ 
120001) . The red channel 3 excess emission could also be 
due to [Fe II]5.34/xm. If H 2 is present in those clouds, 
the IRAC 4.5 /xm band covers the more highly excited 
lines, i.e. the H2 (0,0) S(8)-S(ll) transitions, whereas the 
5.7 and 8 /xm bands cover H 2 (0,0) S(4)-S(7). 

A different type of IRAC emission is mainly charac- 
terized by an excess in channel 2, i.e. it appears as pure 
green emission in Fig. right panel. The faint green 
emission aligned along the southwestern rim is spatially 
coinciding with small clumps of g as seen in the HST 
WFPC2 images ijMorse et al.lll996|r . Judging solely from 
the IRAC band fluxes, it is difficult to determine the na- 
ture of the IR emission in these clumps. Unfortunately, 
none of them are covered by our IRS spectra. Many other 
isolated patches showing channel 2 excess emission do 
not have a counterpart in the narro w band HST WFPC 2 
[O III], [S II], nor the Ha image from lMorse et alJ l|1995|) . 
However, the brightest green patches, one near the cen- 
ter of the ejecta, two blended patches near the north- 
eastern edge of the remnant, and one outside the rem- 
nant perimeter near the northeastern corner of the im- 
age, have optical counterparts visible in archival, broad 
band HST ACS data (proposal 12001, filters F475W, 
F550M, F850LP; see the related Hubble Heritage im- 
age, STScI-2005-30 news release). All counterparts are 
structured, extended sources with color variations. The 
object outside the remnant perimeter is most likely a 
background galaxy, whereas the nature of the objects 
within the perimeter of N132D remains unclear. 

4.2. MIPS continuum emission from shock-heated 
grains 

We interpret the steeply rising continuum in the spec- 
trum of the southeastern rim of N132D as thermal emis- 
sion from dust grains that have been swept up and 
shock-heated by the supernova blast wave. Most of the 
MIPS 24 /xm flux is then due to dust continuum flux, 
since the [O IV] 25.9 /xm line flux is much weaker than the 
continuum in the wide MIPS 24 /xm bandpass (cf. Fig.|SJ- 
The association of the MIPS 24 /xm flux with the blast 
wave along the remnant perimeter and not with the ac- 
tual ejecta (cf. Fig. [21 and 13 confirms that the dust is 
part of the swept-up, ambient interstellar/circumstellar 
material. This is further supported by the appearance of 
the brightest MIPS emission along the southeastern rim, 
where the interstellar gas density is expected to increase 
due to the presence of a molecular cloud (cf. CO contours 
in Fig. [U and Sect. EH}. Generally, the MIPS 24 /xm flux 
intensity follows the X-ray emission very well. Every 
peak in the MIPS 24 /xm surface brightness is also an X- 
ray peak (cf. Fig.EJ), which demonstrates the connection 
between the X-ray emitting plasma and the heated dust 
grains. 
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A two-component modified-blackbody fit to the contin- 
uum yields temperatures of 110 and 58 K (cf. Sect. 13 . 2|) . 
These temperatures can be adopted as average dust 
temperatures if the grains are in thermal equilibrium, 
i.e. balanced heating and cooling rates, and if the mean 
energy per grain heating event is much smaller than 
the thermal energy content of the grain at the equi- 
librium temperatu re (steady vs. stochastic heating, cf. 
IDraine fc Lill200Tl and references therein). Under the 
conditions in the rim of N132D, i.e. an average tem- 
perature T ~ 8. 4 x 10 6 K an d a mean preshock den - 
sity >3cuT 3 (cf. lHwang et~alllT993T: IMorse et al.l fl9961. 
gas- grain collisions probab ly dominate the grain heat- 
ing ijDwek fe Arendd 11992'). If each gas-grain collision 
transfers an energy of the order of kT « 1 x 10~ 16 J, 
i.e. similar to the energy of a 0.6 keV X-ray photon, 
an individual heating event changes the temperature of 
a silicate/graphite grain with radius ^0.04/im by less 
than 1 K. This estimate assumes spherical, homogenous 
grains with a density of 2500 kg m -3 , and a specific heat 
capacity of 185Jkg _1 K _1 for silicate/graphite grains 
at T ~ 100 K. Small grains <0.04/im in size can be 
heated substantially by a single heating event, and thus 
do not have a wel l-defined equilibrium temperature (cf. 
iLi fc Drainell200H) . The fact that we observe two differ- 
ent average grain temperatures may indicate two distinct 
grain size populations and/or emission originating from 
separate physical regions with different heating condi- 
tions (see Sect. 14.31 for a more detailed discussion). 

We can estimat e the dust ma ss Md us t in those two 
components using (Whittet 2003) 

4pF A d 2 



OD\{l dust J 

and adopting a density p = 2500 kg m -3 for sili- 
cate/graphite grains, a distance d = 50kpc for N132D, 
and an average ratio of grain radius a over emissivity 
Qx of 1 X 10 _5 m. We evaluated the flux density F\ and 
the Planck function B\ for the fitted dust temperatures 
Tdust at 23.7 /im. Eq. [5] is an approximation assuming 
spherical dust grains of uniform size, composition, and 
in thermal equilibrium. We derived the dust masses in 
the entire southern rim (outlined by the green color in 
Fig. left panel) by scaling the result from Eq. by 
the MIPS 24 /im flux ratio of the southern rim and the 
IRS slit aperture, which gives Md us t,58K ~ 0.05 M© and 
-Mdust,noK ~ 0.003 M Q . If we further scale this result 
to the total MIPS 24 /im flux of the whole remnant, the 
total swept-up ISM dust mass in N132D becomes about 
0.13 Mq. These numbers have considerable uncertainty, 
particularly the low temperature component, because of 
the assumptions made in Eq. and the measurement 
uncertainties. 

We can compare the derived dust masses with esti- 
mates of the swept-up gas mass from X-ray data, which 
give an approximate upper limit on the total amount 
of dust. iFavata et ail j|1997|) fitted the BeppoSAX X- 
ray spectrum of N132D with a two-temperature non- 
equilibri um ionization mode l. Using the revised best- 
fit from IFavata et all l)1998f) . the total mass of the X- 
ray emitting material is about 380 -850 M w , assuming a 
filling factor between 0.1 and 0.5. IHughes et alJ l|1998f) 
applied Sedov models to fit ASCA spectra of N132D, 
and derived a swept-up mass of about 650 Mq. Assum- 



ing a total swept-up gas mass of 600 Mq for N132D, the 
amount of swept-up dust is ~ 1.5 M Q , if the fraction of 
dust is 0.25 percent of the gas mass, i.e. about one fourth 
of the Galactic value. This dust mass represents an ap- 
proximate upper limit for the case of no dust destruction. 
Considering the factor-of-ten lower dust mass estimated 
directly from Eq. |3 we can state that although some 
portion of the dust in N132D is surviving, as much as 
90 percent might have been destroyed by the blast wave 
shock (cf. Sect. lOjl . 

4.3. Grain processing and destruction in N132D 

Strong supernova shock waves compress, heat, and 
accelerate interstellar material, creating physical condi- 
tions that destroy interstellar dust grains. The grains 
are not instantly following the gas acceleration due to 
their inertia, but are gradually accelerated via gas- 
grain collisions, which cause a gradual gr ain erosion 
through sputtering (cf. iTielens et aflll994j) . Further- 
more, the grain population develops velocity differences, 
and subsequent grain-grain collisions cause shattering 
and vapori zation depending on the relative velocity and 
grain sizes llTielens et al.lll994tlBorkowski fc Dweklll995l 
l.lones et al.lll996D. 

For N132D. lMorse et all l)1996ft have estimated a mean 
shock velocity of the X-ray emitting main blast wave 
of ~800kms _1 , a mean preshock density of ~3cm~ 3 , 
and a postshock t emperature of 8.4 x 10 6 K (cf. also 
iHwang et aflll993|) . which is in good agreement with 
the Sedov model values from IHughes et all |[l998). For 
such a high-velocity shock, the d ominant grain destruc- 
tion process is thermal sputtering iDwek fc Arendtll992t 
ITielens et al"IIT994l iBorkowski fc Dweklll995ft . The life- 
time r of silic ate /graphite grains agains t this process was 
estimate d bv IDra ine fc S alpeterl l)1979f) . and further re- 
fined bv ITielens et alJ (|199^L 



0.79 for graphite 
2.25 for silicate 



x 10 4 yr 



»H 



0.01/mi 



(3) 

at a temperature of 8.4 x 10 6 K. Thus, if the hydrogen 
density is nn = 10cm -3 , silicate/graphite grains with a 
radius a = 0.01 /im are completely eroded after roughly 
1500yr. Likewise, 0.001 /im grains only survive 150yr, 

w hereas . 1 /im g rains would survive 15,00 yr. 

IHughes! l)1987|) . lMorse et al.l ljl996). and H ughes et all 
( 1998) proposed that the supernova of N132D exploded 
most likely within a low-density cavity in order to recon- 
cile the Sedov age with the dynamical age without invok- 
ing large explosion energies. This scenario would allow 
the blast wave to expand unhindered in the beginning 
until it encountered the denser material at the boundary 
of the cavity. If we assume an initial, average expansion 
velocity of ~ 10 kms -1 , the shock reached the cavity 
wall within about 800 years, supposing that the edge 
of the cavity lies close to the current MIPS 24 /im rim 
in the south of N132D and adopting the common cen- 
ter of the oxygen- rich ejecta as the expansion origin (cf. 
IMorse et aljll995|). Since the dy namical age of N132D is 
~2500yr (cf. IMorse et al.lll995f) . we can infer that the 
shock has been interacting with the denser material for 
maybe 1700 yr. Hence, we would expect that most of the 
grains £ 0.01 that were initially present are destroyed 
by now through thermal sputtering if they were continu- 



10 



Tappe, Rho, & Reach 



ously residing in the hot plasma (cf. Eq. [3J . Naturally, as 
the shock wave keeps expanding into new ISM material, 
the grain population is continually replenished. 

In addition to thermal sputtering, shattering and va- 
porization remove some fraction of the large-size grains 
upon grain-grain collision. This is limited by the prob- 
ability of a grain-grain encounter before being brought 
to rest with respect to the gas, which is of the order of 
the dust-t o-gas mass ratio in the absence of bet atron ac- 
celeration (praine fc McKeeHl993tlDrainell2004() . There- 
fore, we conclude that thermal sputtering generally dom- 
inates the grain destruction in the fast shock of N132D, 
but shattering may redistribute the grain size distribu - 
tion appreciably in denser clumps fcf. Uones et alj [l996). 
where the shock speed may be much lower than the es- 
timated mean value of ^800kms _1 . 

We compare our results to calculated model spectra 
fro m dust grains proc essed in fast, nonradiative shocks 
bv IDwek et "all (JT996). Figure ITUI shows the spectrum 
from the southeastern rim of N132D together with our 
modified-blackbody fit components and a s eries of cal- 
culate d shock model spectra adapted from IDwek et al.1 
(1996). Their model assumes silica te /graphite grains in 
a power-law size distribution after Mat his et al.1 l|1977f) 
and grain destruction solely by sputtering without shat- 
tering/vaporization. 

Overall, the scaled model spectra reproduce the ob- 
served spectrum of N132D remarkably well considering 
the model parameters are not fitted to the data. The ap- 
plied scaling factors probably account for the filling fac- 
tor of the emitting material, and were chosen empirically 
to match the flux of N132D. The main difference between 
the observed spectrum and the models is a lack of contin- 
uum emission at 6-9 fim and around 70 /tin, although the 
MIPS 70 /im flux measurement does not provide a very 
stringent constraint due to the background confusion. 
The long-wavelength emission may be due to high density 
material, e.g. a swept-up dense knot, with a shock veloc- 
ity lower than 400 km s~ 1 . The lack of 6-9 /mi flux in the 
models confirms that this emission is due to PAHs and 
not from a population of small, heated silicate/graphite 
grains. Note that the narrow feature near 11.5 /im in the 
calculated spectra is not the 11.3 fim PAH C-H bending 
mode but an IR active mode of crystalline graphite (cf. 
lDrainelli984f) . which would likely be smoothed out i n an 
imperfect polycrystalline sample l)Li fc Draind l2001'). 

Our high-temperature modified-blackbody fit compo- 
nent is very similar to the 1500kms _1 /10cm -3 shock 
model, which justifies that our continuum fit is physi- 
cally meaningful, particularly for the PAH C-C-C emis- 
sion in the 15-20 /jm range (Fig. EJ). Generally, however, 
a single preshock density and/or a single shock velocity 
are apparently not sufficient to explain the continuum 
of the observed spectrum. The medium into which the 
shock wave e xpands is clumpy ( cf. HST WFPC2 images 
of N132D bv lMorse et alJll9 96l. so we expect a combi- 
nation of different densities and shock speeds. 

4.3.1. Comparison with previous observations 

F igure El shows the spectrum of the galaxy NGC 7331 
(cf. iSmith et al.ll2003) in comparison to our spectra for 
N132D. The spectrum of NGC 7331 bears a striking sim- 
ilarity to the spectrum of the molecular cloud/H II- 
region south of N132D, and also shows [Ne III] and 
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Fig. 10. — N132D southeastern rim spectrum with modified- 
blackbody fit components (dashe d lines, cf. Fig.[SJ and calculated 
shock models iDwck ct al. 1996, solid lines). The models are for 
1500/10 (a), 1000/10 (b), 800/10 (c), and 400 kms- 1 / 11 ^ 111-3 (d) 
shock speed/preshock density, and are scaled by factors of 0.13, 
0.14, 0.21, and 0.5, respectively. 

[O IV] /[Fe II] lines. Since it represents an average galac- 
tic ISM mi xture, the resembla nce to our spectra is not 
surprising. ISmith et alJ l)2004|) have fitted the 15-20 /im 
PAH features in this spectrum with two Drude profile 
components at 16.40 (FWHM = 0.34 /im) and 17.12 /tm 
(FWHM = 0.96 /im), and measured a 16.4 + 17.1 PAH 
flux of 1.5 x 10 -7 Wm _2 sr _1 . This fit reproduces the 
overall shape of the emission in this region together with 
the H 2 S(l) component, but it does not take the blended 
features near 17.4 and 17.8 /tm into accoun t, which are 
also o bserved in other sources (cf. Fig. llTJ) . Smith et"al] 
(2004) reported the strong 17.1 /im feature as a possible 
new PAH C-C-C bending mode detection. If we adopt 
Lorentzian components of narrower width, we find the 
same PAH components with similar fluxes and widths as 
seen in our spectra of the molecular cloud/H II- region 
near N132D (cf. Table Hand Fig. [IT]), and no need for 
a broad feature near 17.1 /tm. The 0.2-0.4 /im widths of 
these f eatures are plausible for PAH molecules in the gas 
phase l|Peeters et al.ll2004fl . Our refitted, total PAH 15- 
20 /im flux in NGC 7331 is 1.1 x lO^Wm^sr 1 , which 
is a factor of 1.6 larger than observed in the southeast- 
ern rim of N132D, whereas the PAH band at 11.3 /tm is 
a bout a factor 30 stronge r in NG C 7331. 

iVan Kerckhoven et"aLl (|200fl hereafter VHP) have 
measured 15-20 /tm PAH plateau strengths as well as the 
PAH bands at 6.2 and 11.3 /tm for a variety of Galactic 
sources from ISO SWS spectra. Table[2]shows a compar- 
iso n of d ifferent PAH band ratios from the ISO results 
by IVHPI together with reanalyzed Spitzer archival data 
and our results for N132D. A comparison of the absolute 
fluxes is hampered by the larger ISO beam and different 
source distances, so we restrict ourselves to band ratios 
using Wm~ 2 sr _1 units instead. Note also that IVHP1 
measured the 16.4 /tm PAH features with a local contin- 
uum, i.e. just the excess portion of the feature on top of 
the 15-20 /tm emission hump, whereas we generally use 
the mostly emission-free flux below 15 and above 22 /im 
to determine the continuum for fitting the 15-20 /tm fea- 
tures. We analyze d Spitz er archival calibration data for 
two sources in the VHP] sample, IRAS 22308+5812 and 
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Fig. 11. — Comparison of continuum subtracted 15— 22 fim 
Spitzer IRS emission spectra. All spectra are normalized to their 
16.4 /an flux, and shifted along the flux axis for illustration pur- 
poses. We show our fit for the molecular cloud/H o-region south of 
N132D (online color version: red; cf. Table |T). IRAS 22308+5812 
(compact Hn, high-resolution mode, R ~ 600) has a unique, broad 
feature n ear 19.6 ^m, which also appea rs in ISO SWS data pre- 
sented by Van Kcrckhoven ct al. (2000), and weak emission long- 
ward of 20 urn. See TabldSltor references. 

NGC 7027, in order to test whether the results are com- 
parable. We find that the measured PAH ratios for these 
two sources match within 30 percent, which is a reason- 
able agreement given the uncertainties due to the differ- 
e nt bea m sizes and continuum choices for the fitting. 

IVHPI interpret the observed relative strengths of the 
broad plateaus and the narrower PAH bands in terms of 
variations in the size distribution, i.e. the larger PAHs 
with a few thousand carbon atoms that give rise to the 
15-20 (im emission probably originated from dense cloud 
cores that formed the massive stars in the H II-regions, 
whereas smaller PAHs with less than a few hundred car- 
bon atoms give rise to strong 6.2-11.3/zm bands. They 
also tentatively attribute the discrete 16.4 /jm band to 
C-C-C bending modes of pendant benzene ring substruc- 
tures in smaller PAH molecules (cf. also iMoutou et al.l 
2000, who propose PAH molecules containing pentago- 
nal rings as the carrier of the 16.4 /mi band). 

Table is sorted to show the trend of rising PAH 
15-20 /tm to 11.3 /im ratio. For a diffuse ISM mix of 
starlight-heated, neutral and ionized PAHs, 15-20 /im 
features are efficiently produced by PAH particles in the 
8-25 A size range (300 < A c < 10 4 carbon atoms). The 
11.3 /im feature can be produced by particles as large as 
about 20 A or N c « 4000, but more efficiently by part i- 
cles with radius a < 10 A or N c S 500 llDraine fc L12 0061. 
A direct, quantitative comparison of the observed PAH 
15-20 /im to 11.3 and 6.2 /im, or the 16.4 /im to 11.3 and 
6.2 /im ratios with model calculations is difficult for a 
number of reasons. Foremost, there is a unique PAH size 
distribution, which differs from source to source. Sec- 



TABLE 2 

PAH EMISSION BAND RATIOS IN DIFFERENT TYPES OF SOURCES. 



Source 


Type 


15 — 20 |im 
11.3 Aini 


16.4 (im 
11.3 (im 


15 — 20 lATa 
6.2 (im 


16.4 |im 
6.2 {im 


LkHa 234 H 2 -ridge 


PDR 


0.2 


0.08 


0.1 


0.04 


NGC 7331 (average) 


Galaxy 


0.3 


0.07 


0.2 


0.05 


NGC 7023 


RN 


0.3 


0.14 


0.1 


0.06 


JNCjOyUzy (edge) 


PN 


0.4 


0.04 


1.2 


0.13 


*Nr i r < 7097 


PN 


U . u 




O 7 


u.uo 


*CD-42° 11721 


YSO 


0.7 


0.10 


0.6 


0.08 


*IRAS 03260+3111 


HAeBe 


1.0 


0.06 


0.8 


0.05 


*IRAS 23133+6050 


cH II 


2.2 




1.5 




*IRAS 22308+5812 


cHn 


2.3 




1.8 




IRAS 22308+5812 


cHn 


2.5 


0.27 


0.8 


0.08 


*IRAS 18317-0757 


cHn 


2.5 




1.2 




*Sh 2-106 IRS4 


Hn 


5.1 




3.9 




N132D SE-rim 


SNR 


7.1+1.0 




> 4 




N132D local backgr. 




0.7+0.2 


0.18 


> 0.5 


> 0.13 



Note. — Sources denote d by an asterisk are ISO SWS mea- 
surements taken from Vh3 and converted to Wm — 2 sr — 1 . The 
other sources are Spitzer archival data proces sed and fitte d in th e 
same manner than N132D (NGC 7331, see ISmith eTaTI 120041): 
LkHa 234, see IMorris et all GOOH) : NGC 7023. see IWerner et all 
i200l) : NGC 7027 and IRAS 22308+5812, Spitzer calibration ob- 
servations). Dashes indicate that a discrete 16.4 /an PAH band 
is not detected in that source. PDR = photodissociation region; 
RN = reflection nebula; PN = planetary nebula; YSO = young stel- 
lar object; HAeBe = Herbig Ae/Be star; (c)Hll = (compact)H II- 
region. 

ondly, individual, discrete PAH features may be due to 
specific substructures in PAH molecules of a certain size 
(cf. interpretation of the 16.4 /tm feature above). Finally, 
the specific local physical conditions affect the charge 
state, the amount of dehydrogenation, and the heating 
process of the PAH molecules. 

Nevertheless, it is illustrative to note some general 
trends when comparing the observed PAH emission ra- 
tios in Table [21 with calculated PAH band radiation effi- 
ciencies as function of size for a diffuse IS M mix of neu- 
tral/ionized PAHs (cf. iDraine fc Lil t2006) . We approxi- 
mate the predicted, total emitted power in the 15-20 /im 
region by scaling the calculated 17.4 /tin band power 
with the total integrated absorption cross section per C 
atom in the 15-20 /mi range divided by the cross sec- 
tion for the 17.4 /tm band. This gives approximately a 
scale factor of 4, using the cro ss sections per C atom 
adopted bv IDraine fc Lil (J2006) . The lowest observed 
15-20/11.3, 15-20/6.2, 16.4/11.3, and 16.4/6.2 /tm ra- 
tios in Table [3 correspond then t o the smalles t PAH s 
considered in the calculations by IDraine 
which contain about 20 carbon atoms. PAHs smaller 
than ab out 20-50 C atoms are photolytically unsta- 
ble (cf. lAllamandola et all 119891: iJochims et~aT] 11994b 
lAllain et all I1996H iLe Page et all l2003[) . Hence, the 
PAH size distribution of the sources near the top of Ta- 
ble|3is probably dominated by PAHs of that size. This is 
generally consistent with earlier estimates of PAH sizes 
derived from the 3.3 and 3.4 /im C-H stretching modes 
(|Allamandola et, alJll989l Table 3), and from profile fit- 
ting of the PAH 3.3-11.3 /im bands with calculated spec- 
tra from a distribution of small PAHs ijVerstraete et alJ 
l200lHPech et alJl2002|) . 

For the supernova remnant N132D, we expect that 
small grains/PAHs are rapidly destroyed by ther- 
mal sputtering (cf. Sect. I4.3|> . i.e. relatively larger 
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grains/PAHs tend to dominate the size distribution, 
which is consistent with the high PAH 15-20 /im to 
11.3 /im ratio and the absence of discrete 15-20 /im fea- 
tures. A PAH 15-20 /im to 11.3 /im ratio of ^7 as 
observed in N132D corresponds to PAH sizes of about 
20 A or 4000 C at oms, if we adopt the calculations by 
iDraine fc Lil (|2006) as a tentative guideline. This num- 
ber is a rough estimate, since electron collisions rather 
than starlight dominate grain heating in N132D, chang- 
ing the amount and frequency of energy deposition into 
grains compared to the model calculations. In addi- 
tion, even large PAHs may be dehydrogenated in the 
hot plasma of N132D, which provides another rationale 
for the weak 11.3 /im feature and the resulting high 15- 
20/11.3/im ratio. As a consequence, the size estimate 
from that ratio is probably an upper limit, since PAHs 
in N132D have most likely a higher degree of dehydro- 
genation than in the diffuse ISM. A lower limit estimate 
of the PAH sizes in N132D from the 15-20/6.2 /im ra- 
tio is about 500 C atoms, which is similar to the sizes 
of the very small grains/large PAHs proposed as carriers 
of the bro ad plateaus underlying the 6-14 ^m PAH fea- 
tures (e.g. lAllamandola et alJll98<t iPeeters et"all 12002 
iRapacioli et al.ll2005l and cf. Sect . I3~2~T1 ) . 

The lifetime of 500-4000 C atom grains, about 0.001- 
0.002 /im in size, is 80-160yr according to Eq. [31 which 
describes thermal sputtering of silicate/graphite grains 
in the forward shock of N132D. We presume that large 
PAHs/amorphous carbon particles behave similarly to 
graphite in terms of stability against thermal sputter- 
ing, which is a reasonable assumption since both con- 
sist mainly of aromatic carbon ring systems. Thus, it is 
plausible for large PAHs to survive thermal sputtering in 
N132D long enough to maintain a detectable population, 
particularly since both PAHs and grains will be replen- 
ished as the shock wave keeps expanding into new ISM 
material. 

5. SUMMARY 

We detected strong MIPS 24 /im emission from swept- 
up, shock-heated dust grains in the young LMC super- 
nova remnant N132D (Sect. I4.2fl . In addition, we re- 



ported the first detection of PAH bands from a SNR, 
with tentative PAH features between 6 and 9 /im, a weak 
11.3 /im feature, and a prominent 15-20 /im emission 
hump. This adds to the growing body of sources show- 
ing 15-20 /im emission, which is interpreted as C-C-C 
bending modes of large PAHs (Sect. l3~2~TJl . 

Overall, our observations show that dust grains sur- 
vive a strong supernova shock, which is consistent with 
our estimates of the dust destruction rate by thermal 
sputtering fSect. H3)l . Our results also support the pro- 
posed correlati on of PAH size and the 15-2 to 11.3 /im 
emission ratio l)Van Kerckhoven et al J 1200(1) . The ratio 
of 15-20 /im to 11.3 /im emission in N132D is notably 
higher compared to other types of sources, which we in- 
terpret as rapid destruction of small PAHs/grains in the 
supernova blast wave via thermal sputtering, and pos- 
sible dehydrogenation of the remaining PAH population 
in the high temperature plasma. We conclude that large 
PAHs/amorphous carbon particles with approximately 
500-4000 C atoms survive the harsh conditions in the 
forward shock of N132D long enough to be detected in 
our spectra (Sect. EHTT|) . Future observations of super- 
nova remnants of various ages located in different envi- 
ronments will probe whether, and what fraction of dust 
grains and PAHs are ultimately able to survive the af- 
termath of strong supernova shock waves. 
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